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卫星 轨道 误差 对 定位 精度 影响 的 摄 动 分 析 方 法 
刘 成 ! 李 芳 
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摘要 : 卫星 导航 系统 定位 精度 受 伪 距 测量 误差 、 大气 时 延误 差 、 卫星 原子 钟 钟 差 及 卫星 
轨道 误差 等 多 方面 因素 综合 影响 。 传 统 通常 采用 基于 精度 因子 与 用 户 等 效 伪 距 误差 的 方法 对 
定位 误差 进行 评估 ， 但 在 其 精度 表征 公式 推导 过 程 中 需 对 测量 方程 组 系数 矩阵 H 及 用 户 等 
效 伪 距 误 差 误 差分 布 做 若干 假设 条 件 ， 因 此 它 实 际 上 是 一 个 近似 评估 公式 。 此 外 ， 各 类 误差 
源 中 的 卫星 轨道 误差 属于 三 维 误差 , 需 经 过 坐标 转换 ,并 利用 经 验 参 数 模型 才能 换算 至 用 户 
等 效 伪 距 误差 ,为 此 ,提出 采用 和 矩阵 摄 动 数学 理论 研究 卫星 轨道 误差 对 定位 方程 组 解 的 影响 ， 
利用 谱 范 数 条 件数 对 方程 组 形态 进行 刻画 。 仿真 结果 表明 , 方法 能 够 直接 反映 卫星 轨道 误差 
对 定位 精度 的 影响 , 无 需 进 行 轨 道 坐标 及 用 户 等 效 伪 距 误差 换算 ， 能 够 更 加 直接 和 准确 地 评 
估 卫 星 轨 道 误 差 对 定位 解 精度 的 影响 。 
关键 词 ， 卫星 轨道 误差 ， 导 航 定 位 ; 用 户 等 效 伪 距 误差 范 数 ， 条 件数 
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1 概述 


全 球 卫 星 导 航 系 统 (GlobalNavigationSatelliteSystem,GNSS) 定 位 误差 o, 通常 由 精度 因子 


(Dilution of Precision, DOP) 及 用 户 等 效 伪 距 误差 〈User Equivalent Range Error, UERE) 
m HAZAJ, BU; 
一 ce, = DOP: Oure ， (1) 
A NTRA REPETIR HA EARRA SA PERR KANRA E ROG 
响 的 (统计 〉 和 ， 包 括 伪 距 测量 误差 、 大 气 时 延误 差 、 卫 星 原子 钟 钟 差 及 卫星 轨道 误差 等 。 
在 各 全 球 卫星 导航 系统 建设 与 发 展 过 程 中 ， 上 述 误差 评估 方法 发 挥 了 重要 的 作用 。 它 不 
是 星座 系统 设计 的 重要 依据 ， 也 是 用 户 进行 星座 优选 及 定位 精度 预测 和 分 析 的 重要 指标 。 


Z 根据 精度 因子 参数 类 型 的 不 同 , 用 户 也 能 够 快速 地 对 各 不 同 坐标 方向 上 的 未 知 量 以 及 接收 机 
< 钟 差 进行 解 算 精度 评估 。 
o 然而 , 出 于 方便 使 用 的 考虑 , 基于 精度 因子 与 用 户 等 效 伪 距 误差 的 精度 表征 方法 在 理论 


上 存在 着 一 些 不 足 。 首先 , 精度 因子 与 用 户 等 效 伪 距 误差 的 精度 表征 公式 是 基于 两 个 假设 条 
件 推导 得 到 的 : CIO 假设 观测 方程 组 系数 矩阵 五 无 随机 分 量 ， 这 样 才能 在 推导 过 程 中 将 系 


ZOBEE 五 移 除 到 期 望 算 子 之 外 ;2) 假设 所 有 的 UERE 均 具 有 相同 的 均 方 差 c,， 且 是 互 不 


相关 的 零 均 值 ， 这 样 才能 将 期 望 算 子 简化 为 cz7 。 但 是 ， 只 有 当 所 有 观测 卫星 都 成 正 交 分 


MKRN, MMSR H 才 无 随机 分 量 ， 而 用 户 等 效 伪 距 误差 的 构成 十 分 复杂 ， 
除 符合 高 斯 白 噪声 的 随机 测量 误差 外 , 还 有 各 类 大 气 时 延 、 硬件 设备 造成 的 偏差 ， 以 及 接收 
端 处 的 多 径 效 应 等 误差 。 因 此 ， 在 实际 定位 中 ， 这 两 个 假设 条 件 都 难以 成 立 。 

其 次 , 在 各 类 定位 误差 源 中 , 伪 距 测量 误差 、 大气 时 延误 差 及 卫星 原子 钟 钟 差 等 属于 一 
维 矢量 , 即 卫 星 至 用 户 视线 方向 上 的 误差 ,因而 可 以 轻易 地 将 其 换算 至 用 户 等 效 伪 距 误差 中 。 
然而 ， 卫 星 轨 道 误差 〈 一 般 认 为 是 卫星 广播 星 历 误差 ) 属于 三 维 方向 上 的 误差 量 ， 难 以 直接 
换算 至 用 户 等 效 伪 距 误差 。 为 此 ,为 评估 卫星 轨道 误差 对 定位 精度 的 影响 ,一般 需 将 卫星 广 
播 星 历 误差 从 地 心地 固 坐 标 系 CEarth-Centered Earth Fixed, ECEF) 转换 至 该 时 刻 的 轨道 坐 


标 系 (RTN 坐标 系 ) 下 ， 得 到 卫星 轨道 位 置 在 径 向 CORO. Bü] CT) 和 法 向 CONO 上 的 偏差 
[8。 在 此 基础 上 ， 认 为 径 向 方向 轨道 误差 对 用 户 定 位 精度 的 影响 最 大 ， 再 利用 相关 经 验 模 
型 与 参数 计算 得 到 卫星 URE 值 ， 并 纳入 用 户 等 效 伪 距 误差 中 进行 综合 考量 中。 对 于 我 国 
北斗 卫星 导航 系统 (BeidouNavigationSatelliteSystem, BDS) 等 由 混合 星座 构成 的 全 球 卫星 导 
航 系统 而 言 ， 还 需要 对 MEO、GEO 及 IGSO 等 不 同 轨道 类 型 的 卫星 采用 不 同 的 经 验 参数 ， 
这 增加 了 精度 分 析 与 评估 过 程 的 困难 和 繁琐 性 。 

为 此 , 本 文 提出 采用 和 抢 阵 摄 动 数学 理论 研究 和 分 析 卫 星 轨 道 误 差 对 定位 观测 方程 组 解 的 
影响 ,并 利用 基于 谱 范 数 的 条 件数 对 测量 方程 组 的 形态 进行 刻画 (实际 上 ,精度 因子 的 数学 
本 质 也 即 为 矩阵 范 数 中 的 Frobenius 范 数 )。 仿真 计算 结果 表明 , 基于 矩阵 摄 动 理论 能 够 直接 
反映 卫星 轨道 误差 对 定位 精度 的 影响 , 无 需 进 行 等 效 伪 距 误差 转换 ， 从 而 能 够 更 加 直接 和 准 
确 地 研究 定位 方程 组 性 能 及 其 解 的 误差 ， 以 及 卫星 轨道 误差 对 定位 解 精度 的 影响 。 


2 观测 方程 与 误差 方程 
全 球 卫星 导航 系统 单 点 定位 解 的 线性 方程 组 矩阵 形式 可 表示 为 
HAX =Ap， Q) 


其 中 ， 五 为 方向 余弦 矩阵 ; 


cosa, cos p, cosy, 1 


_|cosæ, cos, cosy, 1 5 


cosa, cos p, cosy, 1 


H}, cosa, cosp, M cosy (i=1,2, ,n) X973 HH P ittama i MEEKER 


量 方向 余弦 。 
AX 为 近似 解 与 真 解 之 间 的 改正 向 量 : 


AX =[Ax Ay Az —cAt], (4) 


Ap 为 由 用 户 估计 位 置 和 钟 差 估计 值 计 算得 到 的 几何 距离 与 观测 伪 距 之 差 : 


Mp=[Ap, o, Ap,] ©) 


lK 


n=4 时 ，(2) 式 为 正定 方程 。 其 解 为 : 


AX=H Ap ， (6) 
当 n>4 时 ，(2) 式 为 超 定 方程 。 采 用 最 小 二 乘法 获得 其 见 余 解 : 
AX -(H'H) H' Ap (7) 
根据 文 [1-3] 中 对 最 小 二 乘 解 算 值 协 方差 窍 阵 的 推导 , 可 得 用 户 坐 标 位 置 误差 的 协 方差 为 
Turca eodd ek n=4 " 
cov(a ) = CH" H)  H' -E[Ap- Ap ]H(H' H)! n>4 


其 中 , EZJBESSET: o, 为 用 户 位 置 误差 均 方差 ，Ay 为 伪 距 观测 误差 向 量 ; m 为 观测 卫星 


rU XE 
C hInaAIV 


x, ES 
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中 
o 


在 上 述 推导 过 程 中 ， 曾 假设 系统 矩阵 H 无 随机 分 量 ， 因 此 可 被 移 到 期 望 算 子 之 外 。 再 进 一 
步 假设 所 有 的 用 户 等 效 伪 距 误差 均 有 相同 的 方 莽 ,， 且 是 互 不 相关 的 零 均 值 ， 则 (8) 式 中 的 
期 望 算 子 可 简化 为 


E[Ap: Ap ]- 6; , (9) 
其 中 ，7 为 4x4 的 单位 矩阵 。 于 是 ，(8) 式 中 的 两 个 方程 可 简化 为 
co(c,)-0;(H'H)'. (10) 


Kb, (HT H) SH P EAE ACE EH P ERANA AKA TRARA, "E 
与 用 户 至 卫星 的 单位 矢量 所 勾勒 的 立体 体积 成 反比 。 

根据 用 户 所 在 坐标 参考 系 的 不 同 , 精度 因子 可 以 表征 不 同方 向 上 的 位 置 误差 。 当 用 户 位 
于 用 户 等 效 伪 距 误差 坐标 系 时 ， 其 实际 含义 为 


XDOP? 协 方差 阵 | 
YDOP* 
(H H) - ; (1D) 
ZDOP 
协 方差 阵 TDOP? | 
其 对 角 线 元 素 的 平方 根 分 别 对 应 三 维 方向 坐标 位 置 及 接收 机 钟 差 的 误差 放大 因子 精度 因子 。 
由 于 在 精度 因子 的 推导 过 程 中 ， 曾 假设 系数 矩阵 H 无 随机 分 量 ， 所 以 (11) 式 中 和 矩阵 


(H'H) 的 协 方差 阵 也 被 认为 是 零 值 。 由 此 可 以 发 现 ， 精 度 因子 的 定义 与 矩阵 范 数 中 的 


Frobenius 范 数 〈 即 F- 范 数 ， 和 矩阵 中 各 元 素平 方 和 的 平方 根 ) 一 致 ， 其 数学 本 质 即 为 矩阵 F- 
范 数 。 这 也 说 明 , 可 以 利用 矩阵 范 数 摄 动 理 论 进一步 开展 对 测量 方程 组 定位 误差 解 的 深入 研 


zx 


5. 
3 轨道 误差 摄 动 分 析 方 法 


从 以 上 误差 方程 推导 可 知 , 基于 精度 因子 的 分 析 方 法 对 测量 值 的 误差 统计 情况 及 系数 矩 
阵 H 的 随机 分 量 情况 均 做 了 条 件 假设 〈 假 设 测量 值 误差 均 符 合 高 斯 分 布 且 相 互 独立 、 系 统 
和 矩阵 五 无 随机 分 量 )。 而 在 实际 情况 下 ， 卫 星 导 航 定位 的 过 程 很 难 符合 上 述 理 想 条件 。 相 较 
之 下 , 基于 范 数 概 念 的 摄 动 分 析 理 论 可 对 测量 方程 组 的 实际 形态 进行 研究 , 并 对 卫星 轨道 误 
差 的 影响 程度 进行 更 加 直接 和 真实 的 评估 。 

当 考 察 卫星 轨道 误差 对 定位 解 的 影响 时 ， 认 为 对 方程 组 左 端 系 数 和 矩阵 H 带 来 摄 动 误差 
ôH ， 并 给 坐标 及 钟 差 未 知 数 向 量 AX 造成 握 动 6 。 即 有 


(H -- óHY(AX - óX) 2 Ap, (12) 
车 


可 以 证 明 ， 由 于 和 矩阵 H 非 奇异 ， 而 由 卫星 轨 位 误差 等 造成 的 摄 动 误差 oH d. BOB 
H +H 仍然 可 以 保持 它 的 非 奇 异性 0]。 因 此 ， 将 (12) 式 减 去 (2) 式 可 得 : 


X =-(H “HAX + H 6HOX), (13) 
从 而 ， 对 (12) 式 左右 两 端 进行 范 数 表达 可 得 : 


Ue 


A 4^ 


EHRTII 


FR 


FI 
[| J 


jax |< ar ajax] + |a "oer ex (4) 


进一步 地 有 
T 
pup, prepa EE Tg 3 
E E (15) 
jax| ^ 1- Ja lon] - mm " 


这 里 , Inr nr T eos. 即 解 的 相对 误差 为 测量 数据 相对 误差 的 | 五 用 五 e. 


Jp |n | 又 称 为 矩阵 的 条 件数 ， 记 为 


cond(H) -|H ||B `|; (16) 
因此 ，(15) 式 可 记 为 

cond (H )* —- a| 
axl 。 四 m 

ax —541 

tT Jon | sr nz | ties HAST KEAN E UE DUIS 
px]... ion PR. - 

lax] Ial 


这 说 明 当 五 具 有 扰动 OH. 时 , 引起 的 解 的 相对 误差 不 超过 五 的 相对 误差 的 cond( 瑟 ) 倍 
piisi 

从 上 述 分 析 可 以 看 出 ， 当 方程 组 中 的 H 具有 扰动 误差 时 ， 方 程 组 解 的 误差 可 以 由 条 件 
Zt cond(H) 决定 ， 其 作用 可 视 为 误差 传递 的 放大 因子 ， 起 到 类 似 于 精度 因子 的 作用 。 条 件 
数 与 所 取 的 矩阵 范 数 有 关 ， 数 学 上 一 般 采 用 谱 范 数 〈 即 2- 范 数 ) 进行 表征 : 


ap [Am (H"H 
cond(H), - |H|,|R |, = RUIT ; (19) 


ieit, cond (H) 又 称 为 矩阵 的 谱 条 件数 DL1244。 


条 件数 cond (H) 同时 刻画 了 方程 组 的 形态 ， 反 映 了 方程 组 解 对 误差 的 敏感 程度 。 当 


cond ( H ) 相对 较 大 时 ， 称 测量 方程 组 为 病态 方程 ， 当 cozd( 召 ) 相对 较 小 时 ， 则 可 称 为 良 


4 仿真 实验 与 分 析 
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为 分 析 和 验证 基于 范 数 的 摄 动 分 析 方 法 的 可 行 性 , 采用 了 实测 与 仿真 相 结合 的 方法 。 首 
先 ， 利 用 全 球 卫星 导航 系统 接收 机 在 某 一 已 知 坐标 的 固定 位 置 测量 了 一 段 历 元 时 长 的 全 球 定 
位 系统 卫星 轨道 坐标 数据 , 并 且 为 了 便于 分 离 和 控制 误差 , 仿真 生成 了 对 应 的 全 球 定位 系统 
伪 距 数据 ， 其 中 含有 lo=3m 的 伪 距 误差。 其 次 ， 为 考察 卫星 轨道 误差 对 定位 精度 的 影响 ， 
在 卫星 位 置 的 X. 了 .Z 三 维 坐 标 方向 上 又 分 别 生 成 和 添加 了 1o=3m 的 轨道 坐标 均 方 误差 ( 因 
此 轨道 的 三 维 误差 为 lc=5.196m)。 通 过 计算 添加 轨道 误差 前 后 的 定位 结果 ， bumps 
道 误 差 造成 的 实际 解 算 误 差 ,， 并 与 基于 范 数 的 摄 动 分 析 方法 估计 得 到 的 误差 进行 对 比 , 5 
如 图 1. 


EMRE OK) 


定位 误差 〈 米 ) 
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图 1 摄 动 分 析 方 法 误差 估计 值 与 实际 解 算 误差 值 对 比 


可 以 看 出 ， 基 于 和 矩阵 范 数 摄 动 分 析 理 论 ， 能 够 不 依赖 误差 分 布 假设 条 件 、 不 通过 轨道 坐 
标 分 解 与 转换 ， 直 接 评 估 和 分 析 卫 星 轨 道 误 差 对 全 球 卫星 导航 系统 定位 解 的 影响 。 并 且 ， 由 
范 数 摄 动 分 析 方法 计算 得 到 的 误差 上 限 均 恰 好 位 于 实际 解 算 误差 之 上 ， 两 者 吻合 情况 良好 
(之 间 差 异 保持 在 数 米 以 内 )， 说 明 范 数 摄 动 分 析 方 法 能 够 准确 、 可 靠 地 对 全 球 卫星 导航 系统 
定位 解 算 误差 情况 进行 刻画 和 评估 。 

5 总 结 


基于 精度 因子 与 用 户 等 效 伪 距 误差 的 精度 表征 方法 在 星座 优选 、 定 位 精度 预测 与 分 析 等 
方面 发 挥 了 重要 作用 。 但 是 , 精度 因子 与 用 户 等 效 伪 距 误 差 的 精度 表征 公式 基于 若干 假设 条 
件 推导 得 到 ， 因 此 严格 意义 上 只 能 作为 近似 评估 手段 ， 在 理论 上 仍 欠 严谨 。 同 时 ,将 所 有 定 
位 误差 源 折算 至 用 户 等 效 伪 距 误差 概念 中 进行 评估 , 也 为 卫星 轨道 误差 对 用 户 定位 精度 影响 
的 评估 工作 带 来 了 一 定 的 困难 ， 增 加 了 评估 过 程 的 繁琐 性 。 

本 文 提 出 基于 范 数 数学 概念 , 采用 矩阵 摄 动 分 析 方 法 对 卫星 轨道 误差 影响 进行 研究 , 提 
供 了 一 种 新 的 理论 和 实用 方法 ， 以 丰 定 和 弥补 传统 以 来 基于 精 度 因子 与 用 户 等 效 伪 距 误差 的 
精度 表征 方法 。 当 已 知 卫星 轨道 误差 先 验 统计 信息 时 , 即 可 利用 方法 对 卫星 轨道 误差 影响 程 
度 进行 预测 和 评估 , 或 在 利用 事后 精密 星 历 计算 得 到 真实 卫星 轨道 误差 后 , 对 卫星 系统 星座 
性 能 进行 更 加 准确 和 便捷 的 分 析 。 该 方法 能 够 避免 对 测量 值 误差 统计 情况 及 观测 方程 组 系数 
和 矩阵 随机 分 量 情 况 的 条 件 假设 , 无 需 进 行 卫星 轨道 误差 的 解析 与 坐标 转换 , 并 且 不 依赖 统计 


意义 ， 能 够 对 单 次 定位 解 的 误差 摄 动 进行 细致 描述 。 因 此 ,能够 更 加 全 面 和 直接 地 评估 与 分 
析 卫 星 轨道 误差 对 用 户 定位 精度 的 影响 。 仿 真实 验 与 分 析 表 明 , 基于 矩阵 范 数 的 摄 动 分 析 方 
法 计算 过 程 简洁 ， 评 估 结 果 准 确 可 靠 ， 具 有 理论 意义 与 应 用 价值 。 
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Studyon Perturbation Analysis Method of the Influence of Satellite Orbit Error on 
Positioning Accuracy 


Liu Cheng!, Li Fang? 
(1. Beijing Institute of Tracking and Telecommunication Technology, Beijing 100094; 


2. National Astronomical Observatories, Chinese Academy of Sciences, Beijing 100190) 


Abstract: GNSS positioning accuracy is affected by pseudo-range measuring error, atmospheric 
time-delay error, satellite atomic clock error and satellite orbit error, et al. Traditionally, the DOP 
and UERE parameters are used to estimate positioning error. However, several assumptions are 
needed to make in the derivation process of precision representation formula, thus it is actually an 
approximate estimate formula. In addition, the satellite orbit error of various error sources belongs 
to three-dimensional error, which requires coordinate transformation and experience parameter 
model to switch to UERE. Therefore, the paper proposes using matrix perturbation 
mathematicaltheory to study the influence of satellite orbit error on positioning accuracy; and 
describes equation form by conditional number using spectral norm. The simulation result 
demonstrates that this method can reflect the influence directly without transformation of orbit 
coordinate and UERE error so that it is more direct and precise to assess the influence of satellite 


orbit error on positioning accuracy. 
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